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Abstract 

We present a model of the gravitational field based on two symmetric tensors. Gravity 
is affected by the new field, but outside matter the predictions of the model coincide exactly 
with general relativity, so all classical tests are satisfied. We find that massive particles do 
not follow a geodesic while massless particles trajectories are null geodesies of an effective 
metric. We study the Cosmological case, where we get an accelerated expansion of the 
universe without dark energy. We also introduce the possibility to explain dark matter with 
5 gravity. 

INTRODUCTION. 

We know that general relativity (GR) works very well at the macroscopic scales pQ. Nevertheles, 
its quantization has proved to be difficult, though. The theory is non-renormalizable, which pre- 
vents its unification with the other forces of nature. Trying to make sense of quantum GR is the 
main physical motivation of string theories [21 [3]. Moreover, recent discoveries in cosmology [I]- [8] 
have revealed that most part of matter is in the form of unknown matter, dark matter (DM) , and 
that the dynamics of the expansion of the Universe is governed by a mysterious component that 
accelerates the expansion, dark energy (DE). Although GR is able to accommodate both DM and 
DE, the interpretation of the dark sector in terms of fundamental theories of elementary particles 
is problematic [9]. Although some candidates exists that could play the role of DM, none have 
been detected yet. Also, an alternative explanation based on the modification of the dynamics for 
small accelerations cannot be ruled out [TQj, ITT] . 

In GR, DE can be explained if a small cosmological constant (A) is present. In early times, 
this constant is irrelevant, but at the later stages of the evolution of the Universe A will dominate 
the expansion, explaining the acceleration. Such small A is very difficult to generate in quantum 
field theory (QFT) models, because A is the vacuum energy, which is usually very large. 

One of the most important mysteries in cosmology and cosmic structure formation is to under- 
stand the nature of dark energy in the context of a fundamental physical theory [T2l [13] - In recent 
years there has been various proposals to explain the observed acceleration of the universe. They 
involve the inclusion of some additional fields in approaches like quintessence, chameleon, vector 
dark energy or massive gravity; The addition of higher order terms in the Einsten-Hilbert action, 
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like f(R) theories and Gauss-Bonnet terms and finally the introduction of extra dimensions for a 
modification of gravity on large scales (See [E]). 

Less widely explored, but interesting possibilities, are the search for non-trivial ultraviolet 
fixed points in gravity (asymptotic safety [15] ) and the notion of induced gravity [H]-[I2]. The 
first possibility uses exact renormalization-group techniques [20]- [23] together with lattice and nu- 
merical techniques such as Lorentzian triangulation analysis [24J. Induced gravity proposes that 
gravitation is a residual force produced by other interactions. 

Recently, in [25, 26J a field theory model explores the emergence of geometry by the sponta- 
neous symmetry breaking of a larger symmetry where the metric is absent. Previous work in this 
direction can be found in [27]-[33]. 

In this paper, we present a model of gravitation that is very similar to classical GR, but could 
make sense at the quantum level. In the construction, we consider two different points. The first 
is that GR is finite on shell at one loop |34j, so renormalization is not necessary at this level. The 
second is a type of gauge theories, 5 gauge theories (DGT), presented in [SUES], which main prop- 
erties are: (a) New kind of fields are created, (pi, from the originals (pi. (b) The classical equations 
of motion of (pi are satisfied in the full quantum theory, (c) The model lives at one loop, (d) The 
action is obtained through the extension of the original gauge symmetry of the model, introducing 
an extra symmetry that we call 5 symmetry, since it is formally obtained as the variation of the 
original symmetry. When we apply this prescription to GR we obtain 5 gravity. Quantization of 
S gravity is discussed in [37] . 

Here, we study the classical effects of 5 gravity at the cosmological level. For this, we assume 
that the universe only has two kind of components, non relativistic matter (DM, baryonic mat- 
ter) and radiation (photons, massless particles), which satisfy a fluid-like equation p = up. We 
will not include the matter dynamics, except by demanding that the energy-momentum tensor 
of the matter fluid is covariantly conserved. This is required in order to respect the symmetries 
of the model. In contrast to [38], where an approximation is discussed, in this work we find the 
exact solution of the equations corresponding to the above suppositions. This solution is used 
to fit the supernova data and we obtain a physical reason for the accelerated expansion of the 
universe within the model: the existence of massless particles. If massless particles were absent, 
the expansion of the Universe would be the same as in GR without a cosmological constant. In 
the Conclusions we speculate on a possible physical mechanism that could stop the accelerated 
expansion and prevent the appearance of a Big Rip. 

We also introduce a possible explanation of DM with 5 gravity. This proposal looks promising, 
but a more detailed work is required to fit the velocity curves of stars around the center of galaxies, 
for example. 

Therefore, in the cosmological regime, we can say that the main properties of this model at 
the classical level are: (a) It agrees with GR, far from the sources and with adequate boundary 
conditions. In particular, the causal structure of 5 gravity in vacuum is the same as in general 
relativity, (b) When we study the evolution of the Universe, it predicts acceleration without a 
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cosmological constant or additional scalar fields. The Universe ends in a Big Rip, similar to the 
scenario considered in [39]. (c) The scale factor agrees with the standard cosmology at early 
times and show acceleration only at later times. Therefore we expect that density perturbations 
should not have large corrections, (d) Apparently, the dark matter could be explain with S gravity. 



It was noticed in [36] that the Hamiltonian of delta models is not bounded from below. Phan- 
toms cosmological models [391 SO] also have this property. Although it is not clear whether this 
problem will subsist or not in a diffeomorphism-invariant model as 5 gravity. Phantom fields are 
used to explain the expansion of the universe. So, even if it could be said that our model works on 
similar grounds, the accelerated expansion of the universe is really produced by a reduced quantity 
of a radiation component in the universe, not by a phantom field. 

It should be remarked that 6 gravity is not a metric model of gravity because massive particles 
do not move on geodesies. Only massless particles move on null geodesies of a linear combination 
of both tensor fields. 



1 5 Gravity Action and Equation of Motion. 

Our modified model is based on a given action S [(f>}, where (pj are generic fields, to which we add 
a new term equal to a 5 variation with respect to the the original fields. We define 8<fii = <pj so 
that we have: 

S[4> t $[ = So[<f>] + K 2 J ^-[flfaz) (1) 

where k 2 is an arbitrary, but small, constant and the index / refers to any kind of indices. This 
new action shows the standard structure which is used to define any modified element or function 
for 5 type models. In particular, let us consider the Einstein-Hilbert Action. This action involves: 

So[g] = [ d A x^g-(-^-R + L M ) (2) 
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S[g] 



j d A x^g- (-j-R + L M + 4 (GT - kT^) g, u + ^ (A^ + A„ iM ) tA (3) 



with k = i k-2 = f|j Lm some matter Lagrangian and where: 

is the energy momentum tensor. The last term in ([3]), depending upon the auxiliary fields A M , 
is needed to impose the condition T^ v v = as an equation of motion in order to implement the 
S symmetry (jSJ) off shell. This term is null in vacuum. So, we have a gravitation model with two 
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fields, the first is just the usual gravitational field g^ u and a second one g^ u = Sg^ v . These fields 
have the following transformations: 



<^9pv = + £o*/;m 

&9p.v = + "I" 9pp£>0,v 9vp^0,p: 4~ 9p,i>,p£,o 

= 6m + ^p£o,m + Vp£o (5) 

where £q and £f are infinitesimal contravariant vectors. The first transformation is the typi- 
cal general coordinate transformations and the second one is obtained using 5 on the first, with 
£i = <5£q. The modified action given by ([3]) was constructed such that is invariant under these 
transformations, so they are the starting point of our model. 

Now, if we vary the fields A M , g^ v and g^ v in (j3J), we obtain respectively the equations of motion: 

T%, = (6) 
= kT^ (7) 

F^^» x D p D x ~g a p = J-p^ (g aP - - X^) - \ (Rr% - W v ) (8) 



with: 



F (nu)(a/3)p\ _ p{(p,j,){af})) v\ + p((pu)(a^)) p,X _ p{{pu)(aP)) g p\ _ p((p\)(af3)) gt tu 

p({ap){jw)) _ 1 (g a ^gP u + g^g 13 ^ — g al3 g^ v ) (9) 

where (fxu) denotes that the fx and v are in a totally symmetric combination. We can see that 
the Einstein's equation do not change. In the application to cosmology, we assume the presence 
of a perfect fluid, so T Q( g = pg a/ 3 + (p + p)U a Up with p and p the pressure and density respectively 
and they are given by (jHJ). U a is the four- velocity of the fluid and fulfills the identity U a U a — — 1. 
Using this energy-momentum tensor we get: 

5 I^l S ali = \g aP (T» a S u ? + T va S^) (10) 
with S aP = ~g afi - X a ^ - X^ a . Using (j7) and (flOjh P is reduced to: 

F {^)W) P X DpDx ~ afj = Qnu^ja, _|_ pjiuyx^ _ gVpjxP\ a p _ I (RVg* — R^g^ — R aV Q^) (11) 

+ 1 (Rrxr + R va Xf - R^X v . a - R ua X». a ) 
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Since is a auxiliary field to implement that T^ u is conserved, we can fix it. So, from now 
we will use the gauge A M = 0. This gauge preserves general coordinate transformations, but fixes 
completely the extra symmetry such that £i M = 0. Finally, the g^ equation is: 



F (^)W)p* DpDx ~ ga p + * _ Rl urg _ jpvffi = q (12) 

With this equation we can completely determine g^ u if we know g^ u . Because, outside the 
sources (T^ u = 0), g a/ 3 oc g a/ 3 is a solution to ffTTj) and we have g a p = rja/3 in the vacuum, our 
boundary conditions will be g a/ 3 = r] a p and g a p = r\ a ^ far from the sources. 

In order to explore the phenomenology of the model, we need to know the equation of motion 
of a test particle. We discuss this in the next section. 



2 5 Test Particle Action. 

To find the coupling of a test particle to a gravitational background field, we follow the prescription 
to construct 5 models ([1]). We know that, in the standard case, the test particle action is: 

Since in our model the modified action is obtained according to ([1]), where in this case <pi = g^, 
the new test particle action is: 

dt = (14) 

where g^ v = g^ u + ^g^u- But (|T3l) and ffT4l are useless for massless particles. To solve this 
problem, we will start from: 

S [x,g,v] = ~ J dt (vm 2 - v^g^x^x") (15) 
where v is a Lagrange multiplier. This action is invariant under reparametrizations: 

x»(t) = a?(t) 
v(t)di = v(t)dt 

i = t-e{t) (16) 

From ffT5l) . we can obtain the equation of motion for v: 
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v = 17 

m 

If we substitute ( ITT]) in ( !T5|) . we recover (IT51) . In other words, (|T5|) is a good action that include 
the massless case. If we substitute (1X51) in ([T]) now, we obtain: 



S[x, g, g, v, v] = ~ J dt (vm 2 - v 1 (g^ v + K 2 g^) x»x v + k 2 v (m 2 + v 2 g fiv x* x x 1/ )) (18) 

This action is invariant under the reparametrization transformations (fT6j) plus v(t)dt = v(t)dt. 
So, f|T8|) is the action that we need to generalize (fT4j) . Two Lagrange multiplier are unnecessary, 
so we will eliminate one of them. The equation of motion for v is: 

_ _ m 2 + v~ 2 (g^ v + K 2 g^ v ) x^x u 

2K,2V~ 3 g a /3X a X 13 

If we now replace f|T9|) in (fT8l) . we obtain the 5 Test Particle Action: 

if ( Tfl^V \ 

S[x, g, g, v) = - J dt \ 3m 2 v - v~ (g^ + K 2 g^ v ) x^x v + - (m 2 + v~ 2 (g^ + K 2 g fiv ) x^x u ) j (20) 

The equation of motion for v is still given by (jT7|) . If we substitute it in (120]) . we obtain (TH]) . 
So, (1201 is a good modified action to represent the trajectory of a particle in the presence of a 
gravitational field, given by g and g, for the massive and massless case. Evaluating m = in (TT5l) 
and (l20l) . they respectively are: 

s!T =0) [x,g,v] = - l -f dtv- x g^x v (21) 



S {m =^ [x, g 1 g,v] = -\J dtv~% u x^ (22) 

with = g^y + n 2 g^ u . The equation of motion for v implies that, in the usual and modified 
case, a massless particle will move in a null-geodesic. In the usual case we have g^x^x 11 = 0, but 
in our model the null-geodesic is g flu x tJ, x 1 ' = 0. 

On the other hand, to describe the massive case in the usual and modified models, we can 
use (fT3j) and (ITlj) respectively. We know that (fT3j) implies that a massive particle will move in a 
typical geodesic with g fiU x fl x l/ = —1, after choosing t equal to the proper time. On the other side, 
the equation of motion of the modified action is more complicated, however we can use the same 
gauge, g^xf = -1. 

To conclude, in the usual case we have: 
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x M + T^aPiP = (23) 

with: 
where: 

g^x^x" = if: m = (24) 
S^i^x" = -1 if : m > (25) 

On the other hand, the equations of motion for a test particle in the modified case are: 

x" + T» ap x a ±P = if: m = (26) 
M>x y + G»pX a xP + N£ Me x a xPxix e = if: m > (27) 

with: 

r % = 2 [ g_1 ] M7 few 3 + g ^>° ~ 

where: 



g^x" = if: m = (28) 
^a^i" = -1 if : m > (29) 

This means that, in our model, massive particles do not move on geodesies of a four-dimensional 
metric. Only massless particles move on a null geodesic of g^. So, 5 gravity is not a metric theory. 



3 Distances and time intervals. 

In this section, we define the measurement of time and distances in the model. 

The geodesic equation, (123]) . preserves the proper time of the particle along the trajectory. Eq. 
( p7|) satisfies the same property: Along the trajectory g f j /V x^x u is constant. So, we define proper 
time using the original metric g^ u : 
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dr = yj-g^dx^dx" = y/-g 00 dx 



(30) 



Following |41j . we consider the motion of light rays along infinitesimally near trajectories, using 
281) and ( 1301) . to get the three-dimensional metric: 



dl 2 = r )ijdx l dx : > (31) 

goo ( BioSjo 
goo V goo 

Therefore, we measure proper time using the metric g^, but the space geometry is determined 
by both tensor fields, g^ v and g^ u . These considerations are fundamental to explain the expansion 
of the universe with 5 gravity. In the next section, we will see this in detail. 



4 FRW and Photon Trajectory. 

To describe the supernova data, we must use the FRW metric. When a photon emitted from the 
supernova travels to Earth, the universe is expanding. This means that the photon is affected by 
the cosmo logical Doppler effect. So the metric g is: 
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R 2 {t)r 2 










\ 



(32) 



R 2 (t)r 2 sin 2 (6) / 



where R(t) is the scale factor that depends on the time parameter t. Assuming an isotropic 
and homogeneous universe, we can use the following ansatz for g: 



( -c 2 A(t) 



9»v 
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B{t) 








B(t)r 2 










\ 



(33) 



J B(t)r 2 sin 2 (6») / 



So, g only has two independent functions, A(t) and B(t), that depend of the time parameter 
t just like R(t). We know that the Einstein's equations lead to: 



kit) 



3R(t) 



(Pi(t)+ Pi (t)) 



(34) 
(35) 



with f(t) = § (t). But now we have additional equations arising from those of g^ v given by 
(TTTj) . These equations are: 

i|fi W -|^5 W -||S( t )-#( t )i W = (36) 

^ (l£) ~ 5py^ (t) ^ 2 5^y^ (() _ -RC*)^*)^^^) - -R(*)-R(t)>l(t) - 2 J R 2 (*).4(*) = (37) 

To solve the system (I34H37I) . we need equations of state which relate piif) and Pi(t), for which we 
take Pi(t) = UiPi(t). Since we wish to explain DE with 5 gravity, we will assume that the universe 
only have non relativistic matter (cold dark matter, baryonic matter) and radiation (photons, 
massless particles). So, we will require two equations of state. For non relativistic matter we use 
PM(t) = and for radiation pu{t) = |p#(£), where we have assumed that their interaction is null. 
Replacing in fl34|) - fl37|) and solving them, we find the exact solution: 



p(t) = p M (t) + p R (t) 

H*n R x(t) + i 

kC^ X\t) 

p(t) = ~p R (t) 



(3f 



3' 

(39) 



HqQ r 



3kC 4 X 4 (t) 

V / ^)TT(X(t)-2) = 3 -^^t-2 (40) 

A{t) = -- ^ 2 (41) 

W n 2 d X(t)(3X(t)+4) 2 



Bit) ~ "N^rf 3X(t)+4 (42) 



X(t) = ^ (43) 

Where R eq is the scale factor at matter-radiation equality, that is pM(t e q) — PR(t eq ), ^is an ar- 
bitrary parameter, C = ^j, and £l R and VLm are the radiation and non relativistic matter density 
in the present respectively, with VLm = 1 — Q R . We know that Q R <C 1, so Qm ~ 1 and C< 1. 
By definition X(t) ^> 1 describes the non relativistic era and <C 1 describes the radiation era. 

Since we have the cosmological solution of the 5 gravity Action now, we can analyze the 
trajectory of a supernova photon when it is traveling to Earth. For this we use ( 128|) in a radial 
trajectory from T\ to r — 0. So, we have: 



- (1 + K 2 A{t))c 2 dt 2 + {R 2 {t) + K 2 B(t))dr 2 = 
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In the usual case, we have that cdt = —R(t)dr. In the 5 gravity case, we define the modified 
scale factor: 



= jm±«M> (44 ) 

1 + K 2 A(t) 



such that cdt = —R(t)dr now. With this definition, we obtain that: 

n=c Lw) (45) 

If a second wave crest is emitted at t — t\ + At\ from r = r 1; it will reach r = at t = i + At , 

so: 

f t0+At0 dt . , 

n = c / (46) 

Jti+Ati ti\l) 

Therefore, for Ati, At small, which is appropriate for light waves, we get: 

At R(t Q ) 



Ah Rfa) 

Since t is the proper time according to (!30|) . we have that 



(47) 



P = 5^ (48 ) 

where u is the light frequency detected at r = 0, corresponding to a source emission at 
frequency u%. So, the redshift is now: 

We see that R (t) replaces the usual scale factor R(t) in the calculation of z. This means that 
we need to redefine the luminosity distance too. For this, let us consider a mirror of radius b that 
is receiving light from our distant source at r\. The photons that reach the mirror are inside a 
cone of half-angle e with origin at the source. 

Let us compute e. The path of the light rays is given by r{p) = pn + ri, where p > is a 
parameter and n is the direction of the light ray. Since the mirror is in r = 0, then p — T\ and 
h = — f i + e, where e is the angle between — ri and n at the source, forming a cone. The proper 



10 



distance is determined by the 3- dimensional metric f[3"Tj) . so we get b = R(t )ri€. Then, the solid 
angle of the cone is: 



/»27r /*e 

= / dcj) sin(0)d0 = 2tt(1 - cos(e)) 



7re - 



;he mirr< 

of all isotropically emitted photons that reach the mirror is: 



where A = nb 2 is the proper area of the mirror. This means that e = — r. So, the fraction 

riR{t ) 



A 

4nr 2 R 2 (t ) 

We know that the apparent luminosity, I, is the received power per unit mirror area. Power is 
energy per unit time, so the received power is P = where hu is the energy corresponding 

to the received photon, and the total emitted power by the source is L = where hv\ is the 
energy corresponding to the emitted photon. Therefore, we have that: 



p = m Lf 

R*(t ) 



R 2 (t )4irr 2 R 2 (t 



where we have used that xt 2 - = — = ^ttt- On the other hand, we know that, in a Euclidean 
space, the luminosity decreases with distance di according to I = t~tf- Therefore, using (14"5"|) . the 
luminosity distance is: 



# 2 (t ) 
— ri 



Moreover, we can use (|40p to change the t variable for Y = CX{t) = (the scale factor 

normalized to one in the present), and then define Y = -|j^y- Using (jH]) and (JI2"]) in (JHJ), we 
obtain: 
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Y\Y Cl]=Y / (3 y + 4 ^) (±CY(3Y + 4C)-lVYTC) 
1 ' ' J V AC (Y(3Y + AC) 2 - IVyTC) 1 } 

We notice that Y ~ Y in the radiation era, that is F <C C, so the universe evolves normally 
in the beginning of the universe, without differences with GR. But when Y ^> C, Y ^> Y and we 

2 

will have a Big Rip, i.e. Y = oo, when the denominator is null. That is Y ~ 5 . We will give 
more detail for this when we will study the supernova data. Now, with all our definitions, the 
luminosity distance is reduced to: 



C Y 2 [C, 



Y'dY' 



H ^Y[Y,C,l}J Y YlY'.C^y/Y^+C 



(52) 



with Y [C, I] = Y[l,C,l}. This means that the distances will be different now. 
In GR, we have: 



c f 1 dY' 

L ~ yh~ Jy v^a^ /4 + n M Y> + n R (53) 

where = 1 — — &r is the dark energy density in the present. We will use (|53|) to compare 
both, a universe with dark energy and our modified gravitation model, with the supernova data. 

Finally, we note that ( flOl) gives us the time coordinate. In the new notation, it is: 



2C; 



R 



(VY~TC{Y-2C) + 2C^ (54) 



Therefore, it is possible to obtain a different age of the universe. A different perception of 
the distances implies a different perception of time. All these differences arise a consequence of the 
modified trajectory of photons. 



5 Analysis and Results. 

Before we analyze the data, we will define the parameters to be determined. In GR, di depends 
upon four parameters: Y, H = lOO/i km s _1 Mpc -1 , Qm and Qr according to (|53l) . However, the 
CMB black body spectrum give us the photons density in the present, f2 7 , and if we assume that 

tt R — Q 7 + fl u — ^1 + 3 (|) (]j) 4 ^ 3 ) fi 7 , we obtain h 2 Vt R = 4.15 x 10~ 5 . Therefore, the parameters 

in di can be reduced to three: Y, h and Ii 2 Qm- For the same reasons, in our modified gravity, di 
depends only on three parameters: Y, C and /, as shown in (!52l) . We use H \fVL R = 0.644 km s _1 
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Mpc" 1 



The supernova data gives the apparent magnitude as a function of redshift. For this reason, 
it is useful to use z instead of Y. So, we have: 

In GR: 

m[z,h,h*tl M ] = ^ + 51og 10 ( dL[ ^p c 2 " M] ) (55) 

, r , , - . c(l + z) Mpc s Z" 1 dY' , . 

100 km ^(^ _ h?Sl M - h 2 n R )Y" + h 2 Q M Y' + h 2 Q R 

YW = ^ (57) 

In our modified gravity: 

m[z,Cj] = M + 5\o gl J dL ] Z ' CJ] ) (58) 

\ 10 pc / 

dikC,Z = c(l + z)- / , 59 

Y[Y(z),C,l) = (60) 

where m is the apparent magnitude and M is the absolute magnitude, common to all super- 
nova, so it is constant. To find Y(z), we must solve ( )60|) using (15 ip . numerically. Now we will 
introduce the statistical method to fit the data. 

We interpret errors in data by the variance a in a normally distributed random variable. If we 
are fitting a function y(x) to a set of points (xi,yi) with errors (cr X i, a y i), we must minimize 



.2 



1 * 

(per point) = — 



(l/i - 2/(^)) 2 



2 2 i // \2 2 



'ft - "yi 

In our case, we want to fit the data (zi, m%) with errors (a Z i, a m i) to the model: 



m(z) = M + 5 log 



10 



lOpc 



Therefore, we must minimize: 



\ 2 



(per point) = If -^^g- (.1) 



i=l u mi ^ \ dz \ Zi n ° 
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Now, we can proceed to analyze the supernova data given in [13] . In both cases, cIl is given by 
an exact expression, but we need to use a numeric method to solve the integral and fit the data to 
determinate the optimum values for the parameters that represent the m v/s z of the supernova 
data. Minimizing (loTj) . we obtain: 

In GR h = 0.6603 ± 0.0078 and h 2 Vt M = 0.096 ± 0.014 with x 2 (per point) = 1.033. 

In our modified gravity I = 34.42±4.27 and C = (3.61±0.21) x 10~ 4 with x 2 (per point) = 1.029. 

With these values, we can calculate the age of the universe. We know that, in the usual case, 
it is 1.37 x 10 10 years, but that in our model it is given by ( 1341 . Substituting the corresponding 
values for /, C and taking Y = 1, we obtain 1.92 x 10 10 years. Finally, we can calculate when the 
Big-Rip will happen. For this, we need Y when Y = oo. That is Y = 1.71, so teig Rip = 4.3 x 10 10 
years. Therefore, the universe almost has lived half of its life. Phantom fields also produce a 
cosmological model that have this property [3HJ ED]. In the S gravity model we can avoid a Big 
Rip at later time by a mechanism that give masses to all massless particles. Some options are 
quantum effects (which are finite in this model) or massive photons due to superconductivity [2] 
which could happen at very low temperatures, which are natural at a later stages of the expansion 
of the Universe. 



6 Dark Matter. 



The principal objective of this paper is to explain the accelerated expansion of the universe with- 
out dark energy, using a modified gravity model. But it is possible that this model explain dark 
matter too. In this section we present a preliminary discussion of this important topic. 

Far from a source, the gravitational field corresponds to a point-like source and we have a 
Schwarzschild solution for g^. With this, we can get a solution to g^ u using (TTTj) . So, we have: 



9u» 



9^u 
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( 



cr(l 



2GM N 



V 
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where b is a new parameter. We have imposed the boundary condition g^ v ~ rj^ and g 



for r — y oo. On the other hand, if we define h^ v 



equation of motion for a test particle ([2 



9iiv ~ Vin> and 
we obtain that: 



Qiiv ~ Vfivi an d we use the 
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2 + K 2 

Where x is the space vector and is the potential. In GR (j) = — therefore = — for 



r — > oo. So, in our model, we will have, for r — > oo, that: 



GM ( K 2 b 



x ~ 2 



Then, if we define the detected total mass, Af r , such that = -£Mt ; we wi u get: 



Mr = M + Mdm (62) 
M DM = (63) 

where M^m can be interpreted as the dark matter. This means that, with 6 < 0, we could 
explain dark matter. We can calculate the dark matter with photons too. For this, we use a radial 
trajectory in a null geodesic, given by (|28|) . for r — > oo. With this, we obtain: 

2 / , 2M(1 + k 2 (1 - b))\ A , , / , , 2M(1 + - 6))\ , 2 



-c 2 {1 + k 2 i * dtr+\l + K2 + ^ ) dr 2 = (64) 

Comparing with the usual case where: 



2M\ , „ / 2M\ , 
c 2 1 ] tit 2 + ( 1 + )dr 2 = (65) 



We obtain that M T = m(i+*2(i-6)) Therefore: 

J 1 + K2 

Af DM = -gL (66) 

1 + K 2 

We notice that photons and massive particles see different M^m, but since K 2 is very small, 
this difference is hard to detect. This means that (1T41) and f l20l) is a consistent system to describe 
dark matter. These preliminary results provide an excellent motivation to explain dark matter 
with this modified model. 
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CONCLUSIONS. 



We have proposed a modified gravity model with good properties at the quantum level. It is finite 
on shell in the vacuum and only lives at one loop. It incorporates a new field g^ v that transforms 
correctly under general coordinate transformation and exhibits a new symmetry: the 5 symmetry. 
The new action is invariant under these transformations. We call this new gravity model 5 gravity. 
A quantum field theory analysis of 5 gravity has been developed [37J [38] . 

In this paper, we study the classical effects in a cosmological setting. To this end, we require 
to set up the following two issues. First, we need to find the equations for 5 gravity. One of them 
is Einstein's equation, which it gives us g^ u , and the other equation is (JHJ) to solve for g^ u . Second, 
we need the modified test particle action. This action, fl20|) . incorporates the new field g^ u . We 
obtain that a photon, or a massless particle, moves in a null geodesic of g^ = g^ u + K 2 g^ v and 
that a massive particle is governed by the equation of motion fl27|) . With all this basic set up, we 
can study any cosmological phenomenon. 

In [38] it was shown that 5 gravity predicts an accelerated expansion of the universe without a 
cosmological constant or additional scalar fields by using an approximation corresponding to small 
redshift. We find in the present work an exact expression for the cosmological luminosity distance 
and verify that 6 gravity do not require dark energy. With this exact expression, we could also 
study very early phenomenon in the universe, for example inflation and the CMB power spectrum. 
This work is in progress. 

On the other hand, photons move on a null geodesic of g, so we can define a new scale fac- 
tor R(t). If we assume that the universe only has non relativistic matter and radiation, we can 
obtain an exact expression for R(t). It is clear in (lo"Tj) that 1 ^> C ^ is necessary to obtain an 
accelerated expansion of the universe. Therefore a minimal component of radiation explains the 
supernova data without dark energy. In this way, in this model, the accelerated expansion of the 
universe, can be understood as a geometric effect. 

We also calculate the age of the universe. We find that the universe has lived a bit more as in 
GR. This is not a contradiction, but rather a reinterpretation of the observations. This result is a 
consequence of the new equation of motion for the photons. This model ends in a Big Rip and we 
calculate when it will happen. The universe almost has lived half of its life. Even though the Big 
Rip could be seen as as a problem, we observe that other cosmological models share this property 
too [3"5l HO]. Nevertheless, in our case, we have some way outs from the Big Rip. For example, 
the appearance of quantum effects or massive photons at times close to the Big Rip, by effects 
similar to superconductivity [33]. These effects could happen at very low temperatures which are 
common at the later stages of the evolution of the Universe. 

Finally, we introduce a preliminary scenario to explain dark matter with 5 gravity. We find 
that, if we assume a point-like source, the Schwarzschild solution is obtained for g and a similar 
solution for g but with extra parameter, b. If we study the asymptotic region, we find that an 
extra mass is perceived by a test particle, massive or massless. This extra mass is determined 
by b and could be understood as a virtual dark matter. To determine whether or not this is a 
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solution of the dark matter problem we have to study the motion of stars around the center of 
the host galaxy. Such a study is been pursued now but goes beyond the scope of the present paper. 
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